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into  account  are  thermal  conductivity,  boiling  of  water,  transpiration  cooling,  and  attenua¬ 
tion  of  thermal  radiation  in  the  soil.  Calculations  were  carried  out  for  a  range  of  parame¬ 
ters.  We  present  an  analytic  fit  to  the  soil  temperature  profiles  for  a  1  MT  explosion  of  200 
ft  HOB  at  TOA  as  a  function  of  eround  range.  .  •  ( 


14  SUBJECT  TERMS 

Pre-Shock  Thermal  Layer 

Soi!  Heating' 


Thermal  Conduction  In  Soil 
Soil  Temperature  Profiles 


15  NUMBER  OF  PAGES 
18 


16  PRICE  CODE 


I  17  SECURITY  CLASSIFICAT'QW 
QP  pCOOfyr 

UNCLASSIFIED 


NSN  7540-280-5500 


CECUHI I  r  i^LMaSIHCAIION 

OF  THIS  PAGE 

UNCLASSIFIED 


19  SECURITY  CL*S3iriuA’..ON 

OF  ABSTRACT 

UNCLASSIFIED 


20  LIMITATION  Ol-  AsSIRACT 


Standard  Form  298  (Rev  2-89) 

Prescrioed  by  ANSI  Sia  239  i 8 
298-102 


SUMMARY 


Calculations  have  been  carried  out  for  the  pre-shock  soil  heating  expected 
from  a  nuclear  explosion  having  a  1  MT  yield  and  a  HOB  =  2000  ft  for  several 
ground  ranges.  Assumptions  made  in  the  calculation  favor  a  lower  surface 
temperature  and  a  rapid  decrease  with  increasing  ground  range.  Additional 
calculations  simulating  flash  lamp  experiments  suggest  that  crucial  parame¬ 
ters  must  be  adjusted  to  obtain  qualitative  agreement.  More  experimental 
data,  as  well  as  calculations,  are  required.  Calculations  having  a  larger  trans¬ 
mission  factor,  a  larger  geometric  factor,  and  possibly  a  larger  soil  absorption 
coefficient  (smaller  mean-free  path),  would  result  in  higher  peak  surface  tem¬ 
peratures  and  probably  would  correspond  to  more  realistic  conditions. 


iii 


Acoesslon  Por 

NT IS  GRAAI 

DT>  C  TAB 

□ 

Unarm  ouncod 

□ 

Justif l nation _ 

By 

Distribution/ 


Availability  Codas 
Avail  aod/or 
Spools! 


CONVERSION  TABLE 


Conversion  factors  for  U.S.  customary 
to  metric  (SI)  units  of  measurement 


To  Convert  From 

To 

Multiply  By 

angstrom 

meters  (m) 

1 . 000  000  X 

E  -  10 

atmosphere 

kilo  pascal  (kPa) 

1.013  25  X  E+2 

bar 

kilo  pascal  (kPa) 

1 . 000  000  X 

E+2 

barn 

meter2  (m2 ) 

1.000  000  X 

E-28 

British  Thermal  unit  (thermochemical) 

joule  (J) 

1.054  350  X 

E+3 

calorie  (thermochemical) 

joule  (J ) 

4.184  000 

cal  (thermochemical  )icm! 

mega  joule/m2(M.I/m2) 

4.184  000  X 

B-2 

curie 

giga  becqucrel  (GBq)’ 

3.700  000  X 

E+l 

degree  (angle) 

radian  (rad) 

1.745  329  X 

E-2 

degree  Fahrenheit 

degree  kelvin  (K) 

TK=(t°f  +  459.67)/  1.8 

electron  volt 

joule  (J ) 

1.602  19  X  E-19 

erg 

joule  (J) 

1 . 000  000  X 

E-7 

erg/second 

watt  (W) 

1.000  000  X 

E-7 

loot 

meter  (m) 

3.048  000  X 

E-l 

lool-pound-force 

joule  (J ) 

1.355  818 

gallon  (U.S.  liquid] 

meter3  (m3) 

3.785  412  X 

E-3 

inch 

meter  (m) 

2.540  000  X 

E-2 

jerk 

joule  (J) 

1.000  000  X 

E+9 

joule/kilogram  (.I/Kg)  (radiation  dose 
absorbed ) 

Gray  (Gy)** 

1.000  000 

kilotons 

tera  joules 

4.  183 

kip  (1000  Ibf ) 

newton  (N) 

4.448  222  X 

E+3 

kip/inch2  (ksi) 

kilo  pascal  (kPa) 

6.894  757  X 

E+3 

klap 

newton-second/m2  (N-s/m2) 

1.000  000  X 

E+2 

micron 

meter  (m) 

1 . 000  000  X 

E-6 

mil 

meter  (m) 

2.540  000  X 

E-5 

mile  (international) 

meter  (m) 

1.609  344  X 

E+3 

ounce 

kilogram  (kg) 

2.834  952  X 

E-2 

pound-force  (Ibf  avoirdupois) 

newton  (N) 

4.448  222 

pound-force  Inch 

newlon-Ometer  (N.m) 

1. 129  848  X 

E-l 

pound-force/inch 

newton/meter  (N/m) 

1.751  268  X 

E+2 

pound-force/fool2 

kilo  pascal  (kPa) 

4.788  026  X 

E-2 

pound-force/inch2 

kilo  pascal  (kPa ) 

6.894  757 

pound-mass  (ibm  avoirdupois) 

kilogram  (kg) 

4.535  924  X 

E-l 

pound-mass-foot2  (moment  of  inertia) 

kilogram-meter2  (kg.m2) 

4.214  011  X 

E-2 

pound-mass/fool2 

kilogram/meler3  (kg/m3) 

1.601  846  X 

E+l 

rad  (radiation  dose  absorbed) 

Gray  (Gy)** 

1 . 000  000  X 

E-2 

roentgen 

coulomb/kilogram  (C/kg) 

2.579  760  X 

F.-4 

shake 

second  (s) 

1 . 000  000  X 

E-8 

slug 

kilogram  (kg) 

1.459  390  X 

E+l 

lorr  (mm  Hg,  0°C) 

kilo  pascal  (kPa) 

1  333  22  X  E-l 

'The  hecquerel  (Bq)  is  Ihe  SI  unit  of  radioactivity;  Bp  =  I  evenl/s. 
'The  Gray  (Gy)  is  Ihe  SI  unit  of  absorbed  radiation. 
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SECTION  1 


INTRODUCTION 


Soil  surface  layers  exposed  to  a  thermal 
flux  of  -102  cal/(cm;s)  are  expected  to  be 
heated  quickly  to  temperatures  in  the  range 
of  several  times  103°K.  Such  an  environ¬ 
ment  is  anticipated  from  nuclear  explosions 
having  low  height-of-burst  (HOB)  at  ground 
ranges  of  interest  to  a  hardened  mobile 
launcher  system.  These  conditions  can  also 
be  created  to  some  degree  in  thermal  layer 
simulators  such  as  the  Thermal  Radiation 
Simulator  (TRS)  or  the  condenser  dis¬ 
charge  simulator  (Flash-Bang). 

In  addition  to  wanting  to  know  the 
properties  of  the  vapor  layer  adjacent  to  the 
ground  produced  by  thermal  radiation, 
modelers  of  the  nuclear  blastwave  desire  to 
incorporate  the  effects  of  dust  pick-up  into 
their  fluid  dynamics  calculations,  including 
the  thermodynamic  consequences.  Conse¬ 
quently,  the  energy  or  temperature  resident 
in  the  first  few  millimeters  of  the  soil 
subject  to  lofting  should  be  specified.  It  is 
expected  that  a  significant  temperature 
gradient  will  be  present  in  this  layer  at  the 
time  of  shock  arrival  so  that  temperature 
profiles  at  various  ground  ranges  are 
needed. 

Several  models  of  the  heating  of  the  soil 
and  adjacent  thermal  layer  have  been 
developed.  Predictions  were  made  for  cases 
of  interest  by  H.  Carpenter,  based  prin¬ 
cipally  on  heat  conduction  into  the  soil.  The 
results  generally  show  surface  tempera¬ 
tures  toward  the  upper  range  of  expected 
results. 


Calculations  were  also  made  using  the 
S-CUBED  THRML  model  showing  much 
lower  surface  temperatures  under  the 
combined  effects  of  thermal  conduction, 
soil  water  boiling/transport,  dust  lofting, 
and  radiation  absorption/scattering  in  the 
dusty  layer.  The  predictions  of  THRML  are 
sensitive  to  the  amount  of  dust  lofted  by  the 
efflux  of  steam  from  the  vaporization  of 
free  or  bound  water  in  the  soil. 

In  order  to  provide  a  detailed  estimate  of 
soil  temperature  we  have  developed  a 
simplified  model  (described  below)  to 
evaluate  energy  transport  in  the  soil.  In 
addition  to  the  process  of  thermal  conduc¬ 
tion  included  by  Carpenter,  water  boiling 
and  transpiration  are  nken  into  account.  In 
addition,  the  thermal  radiation  incident  on 
the  soil  surface  is  deposited  exponentially 
in  the  soil.  While  this  model  is  simpler  than 
THRML.  it  contains  basically  the  same 
physical  processes  in  the  soil;  the  effects 
of  the  thermal  layer  above  the  soil  are 
replaced  by  an  empirically  attenuated  ther¬ 
mal  source  incident  on  the  soil  surface. 

In  Section  2  of  this  report  an  outline  of  the 
formulation  is  presented.  Calculations 
covering  a  range  of  parameters  are  re¬ 
ported  in  Section  3.  In  addition,  this  section 
also  presents  an  analytical  fit  to  the 
calculated  soil  temperature  pro-files  at 
shock  arrival  for  a  yield  of  1  MT  at  a 
height-of-burst  of  2000  ft  as  a  function  of 
ground  range.  Conclusions  and  a  summary 
are  contained  in  Section  4. 
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SECTION  2 


FORMULATION  OF  A  SOIL  HEATING  MODEL 


The  current  model  is  similar  in  concept  to 
the  Analytic  Heat  Transfer  Model  of  H. 
Carpenter  in  that  emphasis  is  placed  on 
processes  occurring  within  the  soil.  In 
addition,  the  model  contains  empirically- 
based  descriptions  of  the  thermal  radiation 
from  the  source  and  processes  taking  place 
within  the  thermal  layer  itself.  In  the 
following  sections  we  outline  the  treatment 
of  the  model  for  these  several  processes. 

2.1  SOIL  ENERGY  TRANSPORT. 

The  calculation  of  temperature  change  in 
the  soil  accounts  for  three  processes; 
thermal  radiation  heating,  thermal  conduc¬ 
tion.  and  water  vapor  advection  (consisting 
of  the  evaporation  and  heating  of  water). 
To  describe  the  evolution  of  this  two-phase 
system  in  generality,  it  would  be  necessary 
to  solve  the  mass,  momentum,  and  energy 
equations  for  each  of  the  interacting  (water, 
water  vapor,  and  soil  matrix)  phases.  We 
have  simplified  this  complicated  formula¬ 
tion  by  making  several  reasonable  approxi¬ 
mations: 

•The  water  vapor  expands  at  constant 
atmospheric  pressure;  it  is  not  neces¬ 
sary  to  solve  a  momentum  equation  for 
the  water  vapor  and  the  mass  equation 
can  be  employed  to  evaluate  the  vapor 
streaming  speed. 

•The  temperatures  of  the  water  and  soil 
matrix  equilibrate  rapidly:  it  is  not 
necessary  to  solve  jeparate  energy 
equations  for  the  water  and  the  matrix. 

•The  soil  matrix  is  not  disrupted  by 
thermal  expansion,  by  the  transport  of 
small  dust  particles  or  by  pressure 


build-up  of  the  water  vapor:  it  is  not 
necessary  to  solve  equations  for  the 
motion  or  density  change  of  the  soil 
matrix  since  dust  lofting  is  neglected. 

The  result  of  the  above  approximations  is 
that  a  single  equation  for  the  total  enthalpy 
change  remains.  It  has  the  form: 


bh 

Bt 

where 

BuHw  B  L.  BT  • 

+ - —  k  —  +  5  (1 

Bx  Bz  Bz 

z 

=  depth  into  the  soil  (measured 
from  the  surface). 

t 

=  time  (measured  from  onset  of 
thermal  radiation). 

u 

=  speed  of  water  vapor  (possibly 
including  partially  unvaporized 
water)  due  to  vaporization  and 
heating, 

H 

=  enthalpy  per  unit  volume  of 
liquid  and  vaporized  water  and 
soil  matrix. 

Hw 

=  enthalpy  per  unit  volume  af 
water  vapor. 

T 

=  common  temperature  of  water 
and  soil  matrix. 

K 

=  thermal  conductivity  of  soil  anti 
water. 

S 

=  volumetric  heating  rate  due  to 
absorption  of  thermal  radiation. 

Equation  ( 1 )  has  the  form  of  an  advec- 
tion-diffusion  equation.  The  advected  en¬ 
thalpy.  however,  differs  from  that  of  the 
time  derivative.  A  complication  arises  in 


relating  the  .nthalpy  to  the  temperature  due 
to  the  phase  change  of  water  and  the 
associated  large  heat  of  vaporization. 

We  expect  that  the  effects  of  water  in  the 
soil  will  be  significant  in  reducing  the  soil 
temperature.  During  the  course  of  soii 
heating  almost  all  of  the  soil  water  in  the 
vaporization  zone  will  be  removed  from  the 
matrix  at  temperatures  corresponding  to 
the  soil  surface  temperature.  The  energy- 
expelled  from  the  soil  in  this  manner  can 
be  a  significant  fraction  of  the  energy 
incident  on  tfr  surface.  We  also  expect  that 
the  obliquity  of  the  thermal  radiation  at  the 
larger  ground  ranges,  as  it  affects  the 
transmission  and  source  deposition,  will 
modify  the  range-dependence  of  the 
results. 

In  the  following  subsections  we  describe  in 
greater  detail  the  three  processes  enumer¬ 
ated  above. 

2.1.1  Thermal  Radiation  Heating. 

The  absorption  of  thermal  radiation  by  the 
soil  matrix  is  approximated  by  a  tempera¬ 
ture-independent  absorption  mean-free 
path.  Q  This  coefficient  must  be  deter¬ 
mined  by  transmission  experiments  or 
calculated  from  the  estimated  index  of 
refraction.  For  our  simplified  model,  the 
thermal  flux.  F,  at  a  soil  depth,  z,  is 
approximated  by  the  single  oblique  ray 
solution: 

dF 

F  =  F()e  -  z/ij.1  Q() )  .  u  =  —  (2) 

dz 

where 

Fa  =  flux  incident  on  the  soil  surface 
(discussed  in  Section  2.2). 

p  =  cosine  of  the  average  angle  of 

incidence  ~  — 

R 


where 

h  =  height-of-burst  (HOB) 

R  =  (h2  +  r2)'  2  =  slant  range  from 
fireball  center  to  the  surface 
station  at  ground  range  r 

Such  a  treatment  would  require  refinement 
if  the  firebali  HOB  w'ere  to  become 
appreciably  smaller  thr-  d.e  radius  of  the 
radiating  fireball  at  times  of  inteiest. 

The  energy  deposited  by  thermal  radiation 
is  calculated  in  the  finite  difference  scheme 
so  as  to  conserve  the  incident  energy.  This 
is  accomplished  by  evaluating  radiation 
tiuxes  at  each  of  the  interzone  boundaries. 

2.1.2  Thermal  Conduction. 

The  thermal  conductivity.  K,  o.  typical  soils 
is  large  enough  to  be  a  factor  limiting  the 
rise  in  surface  temperature  and  modifying 
the  profile  of  soil  temperature  at  the  time 
of  shock  arrival.  Since  the  thermal 
conductivity  of  soils  exhibits  a  range  of 
values,  it  is  desirable  to  explore  the 
sensitivity  of  the  soil  response  to  this 
parameter.  In  audition,  it  is  known  that  the 
thermal  conductivity  depends  on  the  water 
fraction  and  temperature  of  the  soil.  While 
we  have  not  found  data  for  soils  on  the 
temperature  dependence,  it  is  known  that 
he  conductivity  of  SiG?  increases  with 
temperature.  We  have  neglected  any 
temperature  dependence  of  K  at  the  present 
time  due  to  lack  of  data. 

Typical  dependences  of  conductivity  on  soil 
type  and  water  fraction.  Xw,  are  discussed 
by  W.  Sellers,  {Physical  Climatology. 
University  of  Chicago  Press.  1965).  Ac¬ 
cording  to  Sellers.  “ . thermal 

conductivity  is  highest  for  a  soil  containing 
abundant  quartz  and  least  for  a  soil  rich  in 
organic  matter.  It  increases  with  increasing 
moisture  content,  since  water  replaces  ait- 
in  the  soil  pore  spaces,  and  is  higher  for 
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frozen  soil  than  for  unfrozen  soil."  Figure 
1.  taken  from  Sellers,  shows  measured 
conductivity  for  four  soils  at  temperatures 
of  4.4°C  (Fairbanks  Sand).  20°C  (Quartz 
Sand),  and  18°C  (Clay  and  Peat). 

The  figure  indicates  the  volume  fraction  of 
soil.  X-i.  in  the  dry  soil.  Curves  cover  the 
range  from  dry  to  saturated  soils  (although 
this  complete  range  is  not  experienced 
under  field  conditions).  It  is  striking  that  the 
sands  show  a  five-fold  or  greater  increase 
in  K  as  water  volume  fraction  increases  to 
a  typical  value  of  8Cr.  K-values  of  the  sands 
and  clay  are  substantially  larger  than  the 
value  used  by  Carpenter.  We  have 
incorporated  the  curve  for  Quartz  Sand  into 


some  of  our  calculation-,.  !r>.  this  case,  the 
value  of  K  decreases  for  that  part  of  the  soil 
from  which  the  water  has  been  vaporized. 

The  thermal  flux,  <J>,  is  evaluated  at  an 
interface  between  soil  cells  by  a  centered 
difference  approximation  to  the  following 
expression 

0  =  -  K—  . 

ft z 

The  resulting  difference  equations  for  Eq. 
(1)  are  solved  explicitly,  requiring  the 
observance  of  a  time-interval  restriction 
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Figure  1 .  Dependence  of  thermal  conductivity.  K.  on  the  volume  fraction  of  water.  Xw. 
for  four  different  soil  types. 


4 


where 

At  =  time-interval 

Az  =  space-interval 

K  =  thermal  diffusivity 

2.1.3  Energy  Advection  by  Water  Vapor. 

Vaporization  of  water  during  soil  heating 
results  in  vapor  escaping  from  the  soil  into 
the  atmosphere:  it  carries  with  it  the  latent 
and  sensible  heat  it  has  acquired  within  the 
soil.  A  comprehensive  calculation  of  this 
problem  requires  evaluation  of  effects  of 
flow  through  a  porous  medium,  tempera¬ 
ture  equilibration  and  the  erosion  and 
lofting  of  soil  particulates.  While  some  of 
these  effects  can  be  treated  in  a  more 
sophisticated  formulation  than  this  one.  we 
did  not  have  the  time  required  to  develop 
and  test  such  a  formulation  for  our  present 
purposes.  We  believe  that  the  predominant 
processes  are  the  ones  included  in  the 
formulation. 

The  approximations  made  in  this  code  are: 

1 .  Water  vapor  streaming  through  the 
soil  is  able  to  reach  complete 
thermal  equilibrium  with  the  soil, 

2.  The  water  vapor  is  able  to  expand 
rapidly  enough  to  reach  ambient 
pressure  ( 1  atmosphere)  throughout 
the  vaporization  zone  and  to  stream 
freely  into  the  atmosphere  above  the 
surface. 

These  approximations  make  unnecessary  a 
solution  of  a  water  vapor  momentum 
equation  or  the  evaluation  of  an  exchange 
term  between  the  water  vapor  and  the  soil 
in  the  energy  equations.  Further,  assuming 
that  thermal  conduction  and  thermal 
radiation  absorption  take  place  primarily  in 
the  soil,  we  solve  a  mass  conservation 
equation  for  advection  of  the  water  vapor 
and  an  advection-diffusion  equation  for  the 
energy. 


After  some  experimentation  with  numeri¬ 
cal  methods,  we  adopted  an  explicit 
scheme  in  which  the  Eulerian  finite-differ¬ 
ence  mesh  attached  to  the  soil  is  dominant. 
Within  the  zone  of  vaporizing  or  vaporized 
water,  the  following  calculations  are 
performed: 

A.  Form  the  vapor  velocity  due  to 
vaporization  and  heating  in  the 
vapor  region  using  mass  conserva¬ 
tion  and  the  assumption  of  atmo¬ 
spheric  pressure.  Calculate  the 
displacement  of  zone  interfaces  at 
the  end  of  the  time-interval  in 
question. 

B.  Form  the  water  vapor  mass  and 
energy  density  in  these  displaced 
zones  (Lagrangian  phase  of  the 
calculation). 

C.  Interpolate  the  water  mass  and 
energy  density  onto  the  original 
Eulerian  zones  by  a  remapping 
technique. 

D.  Recalculate  the  zone  temperature. 

These  steps  complete  the  advection  calcula¬ 
tion  for  a  cycle;  all  quantities  are  now 
associated  with  the  Eulerian  zones  and  are 
ready  for  the  calculation  of  the  next 
time-cycle.  The  advantage  of  the  above 
procedure  is  that  the  computational  mesh 
remains  associated  with  the  soil  so  that  the 
soil  temperature  profile  is  well  resolved. 

2.1.4  Equation-of-State  for  Water  Vapor. 

Using  the  assumption  that  pressure  is 
maintained  at  one  atmosphere,  water  vapor 
density  data  obey  a  perfect  gas  equation  to 
within  a  few  percent: 

P„  =  5.8  A-  1()-4  £'"/<"!• 
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Assuming  that  the  vapor  fills  the  pore  space 
of  the  soil,  we  obtain  an  effective  water 
vapor  density  taking  account  of  the  pore 
volume  fraction  of  the  soil: 

Pi’  =  ( l  —  Xm  —  X0)  p 0 

where  the  soil  is  composed  of  mineral  and 
organic  components  having 

Xm  =  volume  fraction  of  mineral  com¬ 
ponent. 

X„  =  volume  fraction  of  organic  com¬ 
ponent. 

The  specific  heats  (cal/(cmJ  °K))  of  mineral 
and  organic  components  are  taken  to  be 
0.46  and  0.60.  respectively. 

2.2  SOURCE  DESCRIPTION. 


The  thermal  radiation  flux  component 
normal  to  the  ground  has  been  taken  from 
an  analytic  fit  to  calculated  data  of  the 
RECIPE  code.  This  expression  is 


Fr  =  104-r 


100  YS 


where 


R  \  R  -  50  YS 


f(r)  cal / (cm2 s) 
(4) 


YS  =  Y1  3  =  cube-root-scaled  yield  (Y 
in  kT) 


h  =  HOB  (ft) 


R  =  slant  range  (ft) 
f ( t )  =  2t2/(1  +  t4)  and  r  =  10  t/YS 


yields  to  which  it  is  applicable.  It  appears 
to  fit  the  RECIPE  data  for  Y  =  1  MT 
(Carpenter)  at  peak  fluxes,  and  underpre¬ 
dict  the  flux  at  earlier  times.  These  data  can 
be  replaced  (as  in  the  THRML  code)  by  a 
more  accurate  prescription  for  the  thermal 
flux  using  a  table  of  the  RECIPE  calculation 
of  the  specific  source  and  a  geometric 
model  of  the  fireball  thermal  flux. 

The  factor  accounting  for  the  oblique 
incidence  of  the  radiation  (=  h/R)  is  also  an 
approximation  when  the  fireball  comes  into 
contact  with  the  ground.  As  the  fireball 
after  contact  grows  more  rapidly  into  the 
upper  hemisphere,  the  geometric  factor 
becomes  larger  than  that  of  Eq.  (4).  For  the 
case  evaluated  in  the  sample  calculation  (Y 
=  103  kT.  HOB  =  2000  ft)  the  fireball  is 
expected  to  make  surface  contact  during 
the  later  part  of  the  thermal  pulse. 
Improved  geometric  factors  are  available. 

We  modify  the  source  thermal  flux  by  an 
assumed  atmospheric  transmission  factor. 
T,  due  to  air  absorption  or  scattering,  and 
a  surface  albedo,  A.  due  to  diffuse 
reflection.  Values  of  these  quantities 
employed  in  the  calculations  are  discussed 
in  Section  3.  In  addition,  the  transmission 
factor  is  reduced  by  an  empirical  factor  to 
account  for  thermal  layer  absorption  or 
simulator  geometry.  We  also  permit  the 
surface  to  re-radiate  in  the  IR  due  to  the 
surface  temperature.  Ts.  The  resulting  net 
incident  flux  is 


t  =  time  after  explosion  (s) 


F„  =  FR  T(  \  -  A)-  a  T/. 


Since  this  expression  obeys  cube-root  This  expression  is  used  in  Eq.  (2)  of 

scaling,  it  presumably  is  quite  limited  in  the  Subsection  2.1.1 
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SECTION  3 


RESULTS  OF  CALCULATIONS 


3.1  NUMERICAL  CALCULATIONS. 

The  model  described  in  Subsection  2.2  has 
been  applied  to  several  cases  of  interest  in 
the  modeling  of  the  shock  precursor  of  a 
1  MT  nuclear  explosion.  In  order  to 
compare  results  with  those  of  Carpenter  we 
have  evaluated  soil  heating  at  ground 
ranges  of  3000  to  8000  ft.  using  increments 
of  1000  ft.  for  an  assumed  HOB  =  2000  ft. 
As  a  first  case  we  selected  a  value  of  the 
soil  thermal  radiation  mean-free  path.  ff0 
=  .04  cm.  Each  of  the  calculations  was 
carried  out  to  the  time-of-arrival  of  the 
shock  wave  without  precursor:  these  ranged 
from  -0.7  sec  to  -3.0  sec.  Properties  of  the 
soil  at  earlier  times  corresponding  to  the 
arrival  of  a  shock  precursor  were  also 
available. 

Material  properties  of  the  soil  were  selected 
to  correspond  to  “Quartz  Sand”  of  Figure 
1.  initially  containing  a  water  volume 
fraction.  Xw  =  0.07.  The  thermal  radiation 
source  was  calculated  from  the  prescription 
of  Section  2.2.  which  approximates  the 
results  of  the  RECIPE  Code.  The  radiation 
source  was  reduced  to  account  for  effects 
of  surface  albedo  and  transmission  through 
the  atmosphere  and  lofted  dust.  Each  of 
these  factors  was  assumed  to  be  constant 
(although  it  is  expected  that  they  will 
change  during  the  course  of  the  interac¬ 
tion).  The  albedo  (=  0.25)  was  taken  from 
the  work  of  Carpenter.  This  value  is  typical 
of  undisturbed  light  soils.  Thermal  layer 
transmission  was  estimated  from  the 
experimental  data  of  Versteegen.  ( SAIC 
Flash  Lamp  Tests,  in  Proc.  of  Thermal  Layer 
Review  Meeting.  October  18,  1984). 

Making  use  of  his  calorimeter  measure¬ 


ments  for  Yucca  soil,  we  assume  that  about 
0.5  mean-free  paths  of  attenuation  are 
present  in  the  atmosphere  of  the  flash  lamp 
test  cell.  Since  the  radiation  is  substantially 
normally  incident  in  this  case,  a  greater 
attenuation  may  be  expected  in  the  burst 
geometry  at  larger  ground  ranges.  Conse¬ 
quently.  we  employ  a  time-independent 
thermal  layer  transmission  factor.  T=  exp 
(-0.5  R/h),  where  R  =  slant  range,  and  h  = 
HOB.  Clearly,  more  information  is  needed 
to  clarify  the  dependence  of  flux-reducing 
factors  on  system  parameters. 

Incorporating  the  above  parameters,  calcu¬ 
lations  were  carried  to  the  time  of  shock 
arrival  for  a  series  of  ground  ranges.  The 
peak  surface  temperature  achieved  during 
the  calculation,  and  corresponding  time  of 
peak  temperature,  are  shown  as  a  function 
of  ground  range  in  Figure  2.  At  GR  =  3000 
ft  the  peak  temperature  is  still  rising  rapidly 
at  the  end  of  the  calculation  (0.7  s).  while 
at  the  largest  ground  ranges,  the  tempera¬ 
ture  is  falling  slowly.  These  results  show  a 
rapid  decrease  with  ground  range  primarily 
due  to  the  assumed  decrease  of  transmis¬ 
sion  and,  perhaps,  are  representative  of  a 
weak  precursor  case. 

3.2  ANALYTIC  FITS. 

The  peak  temperature  (TMAX  with  units  of 
°K)  can  be  represented  analytically  as  a 
function  of  ground  range  (GR  with  units  of 
ft)  with  excellent  accuracy  by 


T max  =  300 

4-  2400  exp 


(6) 


GR  -  2500 
3140 
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In  Figure  3  we  show  profiles  of  the  soil 
temperature  at  the  times  of  ideal  shock 
arrival  for  several  values  of  ground  range. 
These  profiles  are  characterized  by  a  rapid 
decrease  from  the  peak  in  a  depth  of  <  1 
mm.  At  the  water  boiling  temperature,  they 
show  a  transition  to  a  much  more  slowly- 
decreasing  curve.  For  convenience  in 
incorporating  these  results  into  hydrody¬ 
namic  calculations,  we  have  obtained 
analytic  expressions  approximating  the 
curves  of  Figure  3.  The  fit  consists  of  two 
branches: 


T  = 

where 


T>  37 3°  A 

(1  +F)  T\,axE 


1  +FE4 


(7a) 


T\iax  =  Peak  surface  temperature  as  a 
function  of  GR  given  bv  Eq. 
(6). 

F  =  0.8/(l  +  2  x  KT4  GR). 

E  =  exp  |  -  1 . 1  x  1 05  z  /( 1 000  + 

GR)|. 

GR  =  ground  range  (ft), 

z  =  soil  depth  (cm). 


T  <  373°A'  (7b) 

T=  300  °K  +  13°exp(-(z-z0)  /  D)  . 


where 


z0  =  0.1 7/(1  +2.5  x  1CT4  GR). 


D  =  .05  + .021 
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The  flash  lamp  experiments  of  SAIC  show 
evidence  of  surface  melting  of  the  soil 
samples.  This  is  suggestive  that  the  surface 
has  achieved  at  least  a  temperature  in  the 
range  of  1000-1500  °K.  In  order  to  compare 
with  the  SAIC  data,  we  calculated  several 
cases  in  which  the  obliquity  of  the  thermal 
radiation  was  neglected  in  the  transmission 


function  and  the  source  deposition.  The  soil 
thermal  radiation  absorption  coefficient 
retained  the  value.  So  =  0.04  cm.  used  in 
previous  calculations.  At  GR  =  7000  ft 
(corresponding  to  the  SAIC  simulation), 
the  resulting  surface  flux  is  ~  28  cal/(cm: 
s).  We  also  used  several  values  of  the  soil 
water  volume  fraction,  Xw,  and  obtained 
the  following  peak  surface  temperatures: 


X 

T  (  °K) 

w 

MAX 

0.035 

1 385 

0.07 

1318 

0.14 

1215 

The  peak  temperatures  display  a  weak 
dependence  on  soil  water  in  this  range. 

The  surface  temperatures  obtained  in  all  of 
these  cases  are  comparable  with  the 
minimum  temperature  indicated  by  the 
data.  There  are  several  effects  which,  in 
fact,  should  result  in  higher  calculated 
temperatures. 

1  .The  peak  flux  attributed  by  SAIC  to 
GR  =  7000  ft  [  >  85  cal/(cm:  s)|  is 
substantially  larger  than  that  given  by 
the  source  prescription  (59  cal/cm;  s) 
used  in  the  calculations. 

2.  We  may  have  reduced  the  flux  too 
much  by  assuming  both  an  albedo  and 
a  transmission  factor. 

3.  The  mean-free  path  may  be  signifi¬ 
cantly  smaller  than  the  selected  value. 

Assuming  that  the  peak  surface  flux  ( for 
a  soil  like  the  Yucca  Flats  sample)  should 
be  approximately  the  observed  48  cal/cm2 
s.  we  should  have  increased  the  flux  by 
almost  a  factor  2.  This  would  increase  the 
peak  surface  temperature  to  approximatelv 
2000°K. 
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TEMPERATURE 


X 


DEPTH  INTO  SOIL  (units  of  cm) 


NB:  For  clarity,  each  curve  is  displaced  from 

that  of  the  next  smaller  ground  range  by 
z  =  0.05  cm. 


Figure  3.  Profiles  of  soil  temperature  for  selected  ground  ranges  from  a  1  MT  explosion 
having  HOB  =  2000  ft.  Corresponding  TOA's  are  in  parentheses. 
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